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DUST ATTENUATION OF THE NEBULAR REGIONS OF z - 2 STAR-FORMING GALAXIES: INSIGHT FROM UV, IR, 

AND EMISSION LINES 

S. DE BARROS^’^’^ N. Reddy^’"^, L ShivaeT’^ 

ABSTRACT 

We use a sample of 149 spectroscopically confirmed UV-selected galaxies at z ^ 2 to investigate the relative 
dust attenuation of the stellar continuum and the nebular emission lines. For each galaxy in the sample, at least 
one rest-frame optical emission line (H(y/[N ii] A6583 or [O III] A5007) measurement has been taken from 
the litterature, and 41 galaxies have additional SpitzerlMWS 24 jim observations that are used to infer infrared 
luminosities. We use a spectral energy distribution (SED) fitting code that predicts nebular line strengths when 
fitting the stellar populations of galaxies in our sample, and we perform comparisons between the predictions of 
our models and the observed/derived physical quantities. We find that on average our code is able to reproduce 
all the physical quantities (e.g., UV /3 slopes, infrared luminosities, emission line fiuxes), but we need to 
apply a higher dust correction to the nebular emission compared to the stellar emission for the largest SFR 
(log(SFR/M 0 yr“^) > 1.82, Salpeter IMF). We find a correlation between SFR and the difference in nebular 
and stellar color excesses, which could resolve the discrepant results regarding nebular dust correction at z ^ 2 
from previous results. 

Subject headings: dust, extinction — galaxies: starburst — galaxies: evolution — galaxies: high-redshift 


1. INTRODUCTION 

Nebular emission (i.e. lines and nebular continuum 
from H II regions) is ubiquitous in regions of mas¬ 
sive star-formation, strong or dominant in the opti¬ 
cal spectra of nearby star-forming galaxies, and present 
in numerous types of galaxies. While several spec¬ 
tral models of galaxies including nebular emission exist 


(e-g-, 

Chariot & Lonshettil 120011: iFioc & Rocca-Volmeranse 

1997 

: Anders & Fritze-v. Alvensleben 20031:IZackrisson et al. 

o 

o 

(N 

), the impact of nebular emission on the determina- 


tion of physical parameters of galaxies, in particular at 
high redshift, has bee n largely neglected until recently. 
IZackrisson et all (l2008l) showed that nebular emission can 
significantly affect broadband photometry; the impact be¬ 
coming stronger with increasing redshift, as the equivalent 
width (EW) of emission lines scales with (1+z). The neb¬ 
ular emission contribution to broad-band photometry at high 
redshift can le ad to overestimates of stellar masses (M^) 
and ages te.g.. ISchaemr & de BarrosI 120091: lOno et al.ll 20 ld: 
ISchaerer & de Barr^l201Qh. 

Until recently, the derived evolution of the specific star- 
formation rate (sSFR=star-formation rate/M^) wi th redshift, 
which indicated an sSFR “p lateau” at z > 2 te.g.. lStark et al.l 
IIoTqI: iGonzalez et aP IMoh . was in conflict with the ex¬ 
pected sSFR-z evolution inferred from hydrodynami cal sim¬ 
ulations (e.g., iBouche et al.l 1201^ iDave et all 1201 Ih . which 
suggest that the sSFR should increase with redshift. Account¬ 
ing for nebular emission can possibly reconcile the obser- 
yations with theoretical expectations of the sSFR evolution 
(IStark et al.ll2Q13l: Ide Barros et al.ll2Q14l: iDuncan et akllloT^ 
ISalmon et al.ll2015h . More generally, an accurate and phys¬ 
ically motivated modeling of nebular emission is necessary 
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to obtain more precise constraints on the physical proper¬ 
ties of high-redshift gala xies, the cosmic st ar-formation his¬ 
tory (ISalmon et al.ll2015h . and reionization (iRobertson et al.l 

l2010h. 

While NIR multi -object spectrographs like MOSFIRE 
(iMcLean et al.ll201^ have enabled measurements of a num¬ 
ber of rest -optical emission lines (e.g., H^, Ha [O III]) up 
to z ^ 3.7 (iHolden et al.ll20l4 iKriek et alJl2015h . these lines 
remain inaccessible at z > 3.7, at least until the launch of the 
James Webb Space Telescope. However, it is possible to find 
empirical evidence of strong emission lines at z > 3.7 by ob¬ 
serving a flux increase in a band which cannot be explained 
with pure stell ar emission and where a strong emission line is 
expec t ed te.g.JCharv et al.l 200 l:IShim et al.l 201 l]: IStark et aP 


120131: Ide Bmos et al.l 12014 ISmit et al.l 120141) . Uncertain¬ 
ties remain relatively large because the inferred strength 
of the lines depends on constraints on the stellar contin¬ 
uum and hence on other inferred galaxy properties (e.g., 
age, dust extinction). One way to avoid this issue is to 
model both stellar and nebu l ar emi ssion simultaneously (e.g., 
ISchaerer & de Barr^ l2009l l 2010 l) . This method must rely 
on several assumptions such as the Lyman continuum es¬ 
cape fraction, the metallicity, and the density and temper¬ 
ature of the ionized gas. These parameters remain largely 
unconstrained at z > 3.7, although there is mounting ev- 
idence that ISM physical conditions eyolve with redshift 
(iKewlev et al.ll 201 ^ iNakaiima & Ouchill2Ql4 ISteidel et'H] 
l2Ql4IShaplev et al.ll2Q15h . 

Until now, most studies on the impact of nebular emis¬ 
sion on galaxy physical properties at z > 3.7 have assumed 
that both nebular and stellar emission suffer the same amount 


2010 

. l2(yil IVanzella et alJ IlM Eml iLabbe et alJ 

| 2010 , 

2013 

: IOeschetal.1 l2013allbl 12014 Ide Barros et al.l 

|2014; 

Duncan et al.l 12014 ISalmon et al.l 120151; iSmit et al.l 120141). 


However, evidence for an excess of nebular color excess rel¬ 
ative to the stellar color excess in local starburst galaxies has 
been provided through the comparison of optical depths be¬ 
tween Balmer lines (Ha and H/3) and the stellar continuum 
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(ICalzetti et al.lll994ll997L[200Qh . This difference in color ex¬ 
cess can be explained by the fact that star-formation occurs in 
regions of high gas and dust column density. UV and optical 
stellar continuum are produced by ionizing and non-ionizing 
stars. On the other hand, the nebular emission arises around 
regions of short-lived massive stars which may not leave their 
dusty birthplaces, while less massive stars can diffuse into less 
dusty regions. This could partly explain why nebular emission 
shows a higher color excess than th e stellar continuum emis¬ 
sion for local star-for ming galaxies (ICalzetti et al.lll9^^ . The 
analysis presented in IWild et afl (1201 Ih also shows that the 
nebular to stellar attenuation is related to the sSFR and axial 
ratio. 

Measuring the relative color excess of the stellar contin¬ 
uum and nebular lines at hig h redshift has been difficult until 
recently (iReddv et al.ll2015h due to the lack of measurements 
of both Ha and H/3 lines for large galaxy samples. One way 
to compare the stellar and nebular color excess is to derive the 
dust extinction of the stellar continuum through SED fitting 
and then compare dust corrected SFRs inferred from emission 
lines (e.g., Ha) with another dust-corrected SFR indicator. 
Under the assumption of a constant star-formation history for 
^100 Myr, SFR(UV) and SF R(Ha) should equilibrate to the 
same value (iKennicuttll 19981) and any remaining discrepancy 
between these two SFR tracers may be due to a different dust 
extinction between the regions of nebular emission and the 
stellar continuum, although there are other possibilities (e.g., 
initial mass function variations). The possible effect of the 
star-formation history must be taken into account when using 
SFR(Ha)/SFR(UV) to infer whether there is a difference in 
the color excess of the stellar continuum and nebular regions. 


2009; 


2014 


redshift (e.g.. 

lErb et al.l l2006al: lEorster Schreiber et al. 

Reddv et al.l 

20101: iMancini et all 1201 ll: iHolden et al. 

Shivaei et al. 

120151). the existence of higher nebular 


color excess rela tive t o stellar color exce ss remains unclear: 
lErb et all (l2006ah and iReddy et al.l (1201 Ol) find that assuming 
a similar color excess between the stellar continuum and neb¬ 
ular regions results in the best agreement between UV and 
Ha-based SFRs, whil e several other studies find that extra ex - 
tinction is required (iKashino et al.l 120131: iPrice et all 120141) . 
A possible explanation for these discrepant results found at 
z ^ 2 is that the trend between the nebular versus stellar color 
excess es d epends on SFR, as was suggested in iReddv et^ 
(120101) and lYoshikawa et af] (120101) . or the nebular to stellar 
attenuation can also depend on the sSFR (I Wild et al. 1 120111: 
iPrice et ^120141) . 

Direct measurements of the B aimer decrement (Ha/Hf3 fiux 
ratio) have been lacking for high-redshift galaxies, until re¬ 
cently, and these observations have been possible because of 
recent new instruments as the Wide-Field-Camera 3 of the 
Hubble Space Telescope, the Fiber Multi Obje ct Spectrograph 
on the Subaru Telescope (iKimura et al.ll201Cil) . and the Multi- 
Object Sp ectrometer for Infra- Red Explor ation on the Keck 
Telesc o pe (iMcLean et al.lll998h . At z ^ 1 .S.lDommguez et al.l 
(12013h . IKashino et all (l2013h and IPrice et al.l (120141) provide 
B aimer decrement measurements for relatively large samples 
of galaxies, but via stacking as t he H/3 line is t y picall y un¬ 
detected for individual galaxies. iKashino et al.l (120131) find 
a higher nebular color excess relative to stellar color excess, 
but the ratio between the t wo is lower that the o ne found for 
local star-forming galax ies (ICalzetti et al .1120001) . The results 
from IPrice et aP (120141) . based on stacked data, are consis¬ 
tent with a difference in color excess between nebular and 


stellar emission, but they also examine possible correlations 
between color excess ratio and the other physical parame¬ 
ters (M^, SFR, sSFR). Interestingly, they find that the dif¬ 
ference in the color excess between n ebular and stellar emis- 
sion increases with decreasing sSFR (iPrice et al .11201 4 Fig¬ 
ure 4). They interpret this correlation with a two-component 
dust model, whereby H ll regions are dustier than the diffuse 
ISM outside stellar birth clouds. The continuum light coming 
from galaxies with high sSFR would be almost completely 
dominated by emission from H ll regions, while the emit¬ 
ted rest-optical continuum from low sSFR galaxies would be 
dominated by less massive stars located outside H ll regions. 
Therefore, for high sSFR galaxies, both nebular and stellar 
emission would be equally attenuated and for low sSFR galax¬ 
ies, nebula r emission would be more attenuated that stellar 
continuu m (iPrice et al.ll20l4 Figure 5). An opposite trend is 
found in IReddv et al.l (12015h with the nebular to stellar color 
excess ratio increasing with increasing sSFR, using the largest 
sample including individual B aimer decrement measurements 
at these redshifts. 

We build upon t hese previous s t udies by using the ap¬ 
proach described in ISchaerer et al.l (120131) . which has been 
al ready used t o analy ze a small sample of lensed galaxies 
in ISklias et H] (l2014l) . Namely, we derive the properties of 
a large number of individual galaxies by comparing several 
observables (broad-band photometry, emission line fluxes, in¬ 
frared luminosities) with predictions from stellar population 
synthesis models that include nebular emission. We consider 
a range of star-formation histories in the SED-fitting in order 
to minimize the number of assumptions going into our anal¬ 
ysis. This study has several aims: the first is to test the abil¬ 
ity of our SED fitting method to reproduce observed quanti¬ 
ties, particularly emission line fluxes. Once our SED fitting 
code is tuned to reproduce the observables, we determine and 
compare correlations among physical properties of individual 
galaxies. Our appr oach is complementary to the one used in 
iReddy et^ (120151) . because we are using SED-fitting code to 
examine these issues. 

The paper is structured as follows. The selection procedure, 
spectroscopic and photometric data, and the final sample used 
in this work are described in Sect. [21 The method used to de¬ 
rive physical parameters and to model the nebular emission is 
described in Sect. O The ability of our method to reproduce 
emission line fiuxes is described in Sect. [4] and the implica¬ 
tions regarding the dust extinction toward nebular regions are 
discussed in Sect. [5] Section [6] summarizes our main conclu¬ 
sions. We adopt a A-CDM cosmological model with Hq = 70 
km s ~^ Mpc~^, Um = 0.3 and = 0.7. We assume a Salpeter 
IME (ISaloetedfl9551)^ All ni agnitudes are expressed in the AB 


2. DATA 

2.1. Selection, optical and near-IR data 

All the data used in this paper are taken from the literature. 
In the following, we provide only brief descriptions of the data 
acquisition and reduction. The reader is invited to consult the 
references given below for more details. 

We start with a large UV-selected galaxy sample, 
with phot ometric and spect ro scopic obse r vation s de- 
scribed in ISteidel et alJ (120031) . ISteidel et~n] (120041) . and 
lAdelberger et al ] (l2nnl . The selection is based on UnGlZ 
colors. All the galaxies in our sample are spectroscopically 
confirmed with rest-frame UV observations conducted with 
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Figure 1. Redshift distributi on of the star-forming galaxy sample. All g alax- 
ies have Ha measurements JErb et al.ll2006al : IKulas et 120121120131) . ex- 
cept 5 galaxies at z > 3 for which we have [O iii] A5007 m easurements 
JKulas et ^120121) . [N ii] A6583 measurements are taken from IKulas et all 

the Keck/LRIS spectrograph (lOke et al.lll995h . 

Photometry in the J and Kg band are obtained with the Palo- 
mar/W IRC and Magellan/PANIC instruments (IShaplev et alJ 
120051) . These data provide important constraints on the 
B aimer break, which can be used to estimate the age of the 
stellar population. 

Ou r photometry inc ludes Huhhler^¥C?> F160W observa¬ 
tions (iLaw et al.ll2012h which provide more robust constraints 
on the Baimer and 4000 A breaks. We combine ground- 
based UnGlZ optical photometry with /, Kg, and F160W data 
with Spitzerf[RAC observations in the four available channels 
(3.6/im, 4.5/im, 5.8/im, and 8.0/im). Details r egarding the 
fields o f vie w and photometry ca n be found in iReddv et^ 
(l2006al) and iReddv et al.l (l2012bb . The latter reference pro¬ 
vides a complete description of the data used in our study. 

To more robustly constrain the physical parameters of 
galaxies studied in this paper, we add available near-infrared 
spectroscopic data taken with t he NIRSPEC in strument 
(iMcLean et al.lll99^ presented in lErb et al.l (l2006bh . These 
observations provide Ha fiux measurements for 91 of our 
galaxies. Addition al NIRSPEC spectroscopy is taken from 
IKulas et ^ (l2012h . which provides Ha and [O III] A5007 
fiux measurements for 5 and 5 gal axies respectively. We 
also use spectroscopic observations (IKulas et al.l 120131) per¬ 
formed with the recently commissioned Kec k/MQSFIRE 
near-infrared spectrometer (iMcLean et al .112013) . This allows 
us to add 48 measurements of both Ha and [N ii] A6583 
fiuxes to our study. The main source of fiux measurement 
uncertainties for both NIRSPEC and MOSFIRE instruments 
come from slit losses and are estim ated to be close to ^ 2 
(lErb et al.ll2006al:[Steidel et al.ll2Q14l) . This average slit loss is 
mainly due to seeing conditions and not to the size of galax¬ 
ies: for the stellar mass range probed here, at z ^ 2, half-light 
radius is typically 0'.'2-0'.'3 (iPrice et all 120151) . well within 
the NIRSPEC and MOSFIRE slits used to obtain the data 
(iRrh et al.ll2nn6hl:IKulas et ani2?)ni2(yr^ . Therefore we cor¬ 
rect the observed fiuxe s by a factor 2. The lErb et al.l (l2006bl) 
and IKulas et al.l (120121) samples probe similar range of stellar 
masses and SFRs as stated in th is latter reference. The galax¬ 
ies from the iKulas et ^ (120131) sample are more massive on 
average (IO^^Mq <M*< lO^^M©) but still in the stellar mass 


range probed in lErb et ^ (l2006bl) and iKulas et al.l (120121) . 

Einally, combining photometric and spectroscopic data de¬ 
scribed in this section, we obtain a sample of 149 galaxies 
with extensive broad-band photometric coverage, and Ha and 
[O III] A5007 fiux measurements for 144 and 5 galaxies re¬ 
spectively. [N II] A6583 measurements exist for a subsample 
of 48 galaxies with Ha observations. 

2.2. Mid-IR data 

MIPS 24/im data in the GOODS-North fi eld and four LBG 
fields (Q1549, Q1623, Q1700, and 02343: iRZfdv fc Stei dell 
l2009h are used to measure the rest-frame 8 fim fluxes of galax- 
ies at z ^ 2. T he complete photom etry and extraction method 
is described in iReddv et al.((l201Ql) . The measured 8 jim lumi¬ 
nosities are then converted to the total IR lu minosity (i.e., IR 
lumino sity integrated over 8-1000/im) using IChary & Elb^ 
(1200 ll) templates, and are corrected for the luminosity depen- 
dance of the conversion, as derived from stacked Herschel 
data for a similarly-selected sample of galaxies (Reddy et al. 
(2012a). This latter reference compared L(IR) computed from 
24, 100, 160/im, and 1.4 GHz fiuxes with L(IR) derived from 
the MIPS 24/im data only: L(IR) estimates are consistent with 
each other, within the uncertainties of the measured IR lumi¬ 
nosity. We have 41 galaxies with MIPS 24/im data: 21 with 
detections (S/N > 3) and 20 with upper-limits. The redshift 
distribution of our final sample and subsamples are shown in 
Pigure[I] The constraint provided by the IR luminosity allows 
to compute the bolometric SER (assuming energy balance be¬ 
tween UV and IR), and therefore the dust extinction, in an 
independent way. We provide more details on how we con¬ 
strain galaxy physical parameters with IR data in Section [3l 

To aid our analysis, we define three different subsamples: 
the “MIPS” sample (galaxies with MIPS data), the “MIPS de¬ 
tected” sample (galaxies detected at 24/im with > 3 sigma), 
and the “MIPS upper limit” sample (galaxies undetected), 
with 41, 21, and 20 galaxies respectively. All these galaxies 
have at least Ha measurements, with 12 having [N ii] A6583 
measurements. 

3. SED FITTING 

Physical properties of the galaxies are inferred by fit¬ 
ting model SEDs to the broad-band photometry from U- 
band to 8/im. Our SED fitting code is a modified version 
of th e photometric redshift code HyperZ (iBolzonella et al.l 
I 2 OOOI) . We generate a set of spectral t e mplat es with the 
GAL AXE V code of iBruzual & Charlotl (120031) . for three 
different metallicities (Z=0.004, 0.04, 0.02) and three dif¬ 
ferent star-formation histories (SEH=constant, exponentially 
rising, exponentially decreasing). The timescales used 
for the rising and declining SEHs are respectively = 
[10,30,50,70,100,300,500,700,1000,3000] Myr and rj = 
[10,30,50,70,100,300,500,700,1000,3000, 00 ] Myr. The 
stellar age is defined as the age since the onset of star- 
formation. Unless stated otherwise, we consider a mini¬ 
mum age of 50 Myr corresponding to the typ ical dynami¬ 
cal tim escale expected for z ^ 2 galaxies (e.g., iReddv et^ 
12012H) . 

We assume that the dust att enuation is describ ed by the 
Calzetti dust attenuation curve (ICalzetti et al.ll200Ql) . with Ay 
allowed to vary from 0.0 to 6.0, in steps of 0.05. The inter- 
galac tic medium (I GM) is treated using the Madau prescrip¬ 
tion (lMadaulll995l) . Redshifts are fixed to the spectroscopic 
ones. Eor our entire grid of models, we compute the ^nd 
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Figure 2. Comparison between log SFRuv+ir and log(SFRuv+iR/SFRsED) for three different SFHs (constant, rising, and declining, from left to right). An 
increasing inconsistency between SFRuv+ir and SFRsed appears at SFRuv+ir > 10^-^^ M 0 yr“^ for galaxies detected in MIPS 24/rm. We remind that we assume 
energy balance in our SED fitting procedure. 


the scaling factor of the template; the latter determines the 
SFR and M^. 

We account for the effect of nebula r emission (both con¬ 
tinuum and emission lines) following ISchaerer & de BarrosI 
(l2009l) and ISchaerer & de Barr^ (l201Ql) . The strength of 
nebular emission depends mainly on the number of Ly¬ 
man continuum photons, which is computed from stel¬ 
lar population synthesis models. Relative lin e inten sities 
are taken from lA nd ers & Fritzelv. AlvenslebenI (120031) and 
IStorey & Hummeil (Il995l) . for typical ISM conditions (rie = 
100 cm“^, T = 10"^ K). We assume that the Lyman contin¬ 
uum escape fraction is equal to zero, which means that our 
models produce the maximum theoretical strength for nebu¬ 
lar emission for the adopted ISM physical conditions. Addi¬ 
tionally, we assume that the same dust attenuation curve ap¬ 
plies to the nebular and stellar emission, and that the color 
excess derived from the continuum applies to the line emis¬ 
sion. If, as in the local universe, the nebular lines are subject 
to larger color excess than UV continuum, then our assump¬ 
tion would over-predi ct the emission line fluxes compared to 
the measured fluxes. iReddv et^ (1201 Ql) show that the im¬ 
pact of nebular emission on ages and stellar masses at z ^ 2 
is insigniflcant because most of the longer wavelength bands 
(e.g., IRAC) used in the SED fitting are unaffected by line 
emission at these redshifts. Therefore our assumptions about 
nebular emission modeling, particularly about nebular atten¬ 
uation, does not significantly affect the physical parameters 
derived from the SED modeling. 

We also ignore the Lya line contributions to the broadband 
fluxes because over 122 galaxies with hya flux measurements 
available in our sample, only one third (44) exhibits Lya in 
emission, and more than a half of this third have EW(LyQf)< 
10 A. Therefore we set the Lya flux to zero in our nebular 
emission modeling. 

We used the IR luminosity as an additional constraint to 
the SED fitting by assuming a balance between the energy 
absorbed in the UV/optical (we integrate over 912A to 3/im), 
and re-radiated in the IR. HyperZhuilds a grid of models (Hy¬ 
percube). Eor each model in our grid, we compute the ex¬ 
pected IR luminosity (8-1000/im) and then exclude from the 
grid all the solutions for which L(IR)sed is inconsistent with 


L(IR) 24 ^m within 68% confidence. Eor the non detections, we 
exclude solutions inconsistent with the L(IR) 24 ^m 3 a upper 
limits. The 41 galaxies in the MIPS sample will have the most 
robust determinations of SERs and dust reddening given the 
independent information provided by the IR constraints. We 
have inspected the SED for each object to ensure that the fits 
are still good when assuming the energy balance: there is no 
significant difference between with and without energy bal¬ 
ance (comp ariso n based on x^), except for a small subsample 
(see Section l4Jl) . 

Minimization of over the entire parameter space yields 
the best-fit parameters. Uncertainties are determined from 
the distribution, through the likelihood marginalization for 
each parameter of interest with C oc exp(-x^/2). 

4. RESULTS 

In this section, we first compare t he S ERs derived through 
SED fitting with SERuv+ir (Section [4T]) and determine if th e 
assumed Calzetti reddening curve is appropriate (Section[42]). 
Because our final sample has been obtained after applying 
several criteria (UV-selected, optical emission line measure¬ 
ments), we compare the general properties of our sample with 
previous studies to ensure that it is representative of z ^ 2 star¬ 
forming galaxies (Section [4]3]). We then test the ability of the 
SED models to reproduce observed emission line fluxes for 
the MIPS detected sample (Section 14.4b and explore if there 
is any trend among the physical parameters that could explain 
the discrepant results regarding the difference in color excess 
between nebula r line s and stellar continuum found in previous 
studies (Section [43]). 

4.1. SFRsed versus SFRey+ir 

Eigure [2] shows the comparison of the SERs derived from 
SED fitting with those computed from the sum of the ob- 
scu md (IR) and un obscured (UV) SERs (the latter assume 
the iKennicuttll 19981 relations). We And a good agreement be¬ 
tween SERsed and SERuv+ir for SERuv+ir^ lO^MQyr ^ for 
the rising and constant SEHs: for the MIPS detected sam¬ 
ple, the difference is < 0.25 dex. There is a mismatch be¬ 
tween SERsed and SERuv+ir for 7 galaxies which are the only 
known ultra-luminous galaxies (ULIRGs, L(IR) > IO^^Lq) 
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Figure 3. Relative probability for the SMC extinction curve to be the best-fit 
dust extinction curve versus age (also derived assuming an SMC curve). We 
show the result assuming a constant SFH for the entire sample. The dashed 
lines show the threshold we adopt when determining whether a model is the 
best fit one (pr >0.9). All galaxies above the upper line are considered best 
fit with the SMC curve and all galaxies below the lower line are considered 
best fit with the Calzetti curve. We introduce a slight offset in age to make 
the figure clearer. 


in our sample. For these galaxies, the discrepancy between 
SFRsed and SFRuv+ir is >0.5 dex, and can be as large 
as ^ 1 dex. We further investigate the ULIRGs physi¬ 
cal properties in Appendix O The declining SFH leads to 
lower SFRsed in comparison with SFRuv+ir- This is con¬ 
sistent wit h the fact that declining SFHs can underpredict 
SFRs (e.g.jReddv et al.ll201Ql:IWuvts et ^11201 lb iReddv et HI 

l2012btRice et al.ll2014l) . In the following, we will show and 

discuss the results for the rising SFH mainly, since the con¬ 
stant and rising SFHs lead to similar results (except as indi¬ 
cated otherwise). 

In the remainder, we present results based on SED fitting 
assuming energy balance, but we exclude these 7 ULIRGs 
from our analysis. 


4.2. Attenuation curve 

While we assume a Ca l zetti a ttenuation curve to perform 
SED fitting, iReddv et ^ (1201 Oh show that the IRX-/3 rela¬ 
tion i s more consistent wit h an SMC-like curve (iPrevot et al.l 
11984 iBouchet et al.ll 1985b than the Calzetti curve for galax¬ 
ies younger than 100 Myr. To check if this is also true for 
our sample, w e use the Akaike Information Criterion (AIC; 
lAkaikel[l974l) which allows us to compare different models 
and derive the relative probability for one model to be the best 
model among a given set of models. In our case, we compare 
modeling assuming the Calzetti curve with those that assume 
the SMC curve (where all other SED parameters have been 
fixed in the fitting). AIC is defined as: 

AlC = x>2k ( 1 ) 


where k is the number of free parameters. From the AIC, we 
derive the relative probability pr for the ith model (model as¬ 
suming SMC/Calzetti attenuation curve) to be the best model 
as: 

A, = AIQ-min(AIC) (2) 


exp(-l/2A/) 

i;/exp(l/2A0 


(3) 


We show the results obtained for the entire sample (Figure 


Q1623-BX458 



Figure 4. Observed SEDs (thick black crosses) and best-fit assuming a con¬ 
stant SFH with a Calzetti attenuation curve (in red) and a SMC attenuation 
curve (in blue) for a galaxy best fit (i.e. pr > 0.9) with an SMC curve (top) and 
a galaxy best fit with the Calzetti curve (bottom). Thin red and blue crosses 
show integrated fluxes in each band. The main physical properties are shown 
in blue and red assuming the SMC and the Calzetti curves respectively. Apart 
from the use of SED fitting to discriminate among attenuation curves, it is 
worth noting that the Kg band excess for the galax y Q1623-BX458 can only 
be explained by Ha emission dShivaei et aklDOlST) . 


0. Since AIC is only a statistical test, we define a thresh¬ 
old where we consider a model as the best fit one if the rel¬ 
ative probability is pr> 0.9. While galaxies best fit with the 
SMC curve are found at any age, galaxies of a young age 
(< 10^-^ years) are more likely to be fit with a SMC-like at¬ 
tenuation curve, while galaxies best fit with the Calzetti curve 
are ^10^ years. This relation between age and attenuation 
curve is similar to previous results in literature based on the 
comparison between UV to IR luminosity ratio and 0 slope 
(IReddv et all \20^ IMol 1^123: iSiana et all [200^ IMm 
IW uvts et al.ll2012l:[Shivaei et 70120151). 

We show in Figure [4] two examples: one galaxy best fit 
with the SMC attenuation curve and one best fit with the 
Calzetti curve (pr(SMC) > 0.99 and pr(SMC) < 0.01, respec¬ 
tively). The main difference between the solutions obtained 
assuming different attenuation curves is the UV (3 slope (/a oc 
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Figure 5. Star-formation rate versus stellar mass for three different star-formation histories. Star-formation rates are inferred from SEP fitting assuming en ergy 
balance when MIPS lAfim data are available (see Section [3}. The red and blue lines show the sta r formation main se quence from IDaddi et al.l (120071) and 
IReddv et'di (I2012br) . respectively, at z ~ 2. The blue dashed line shows a linear fit to the data from IReddv et al.l (I2012br) . We show linear fits to our data for 
three different samples: MIPS detected (dashed lines), MIPS data (dash-dot lines) and the entire sample (dotted lines). Typical error bars are shown in the top 
left corner for each SFH. 

our sample. If we consider probabilities at face values, ^ 20% 
of the galaxies are best fit with an SMC curve and ^ 80% 
are best fit with a Calzetti curve. Because the majority of 
galaxies in our sample appears to be better described with the 
Calzetti curve, we assume this attenuation curve in the rest of 
Section[4l 



Figure 6 . UV ^ slope distributions for the IReddv et all (120151) sample, our 
entire sample, and the MIPS sample. The /3 slopes are derived from the 
observed SEDs. 

A^). The UV /3 slope depends on the age and dust attenua¬ 
tion of the galaxy ( and the SFH, iLeitherer & Heckm^l 19951: 
iMeurer et all Il999h . In the two examples shown, Baimer 
breaks are constrained by the J-IZ color but in each case, only 
one attenuation curve is able to reproduce both the B aimer 
breaks and the UV [3 slopes. Best fit assuming the Calzetti 
curve leads to /3 slope bluer than the slope derived with the 
SMC curve. For young galaxies (< lOOMyr), the Calzetti 
curve leads to (3 slopes too blue compared to the observed 
ones (Figure [4] top), while for old galaxies, the SMC curve 
leads to /3 slopes too red (Figure [4] bottom). Despite the well 
known degeneracy between dust attenuation and age deriva¬ 
tion, the color is enough to lift the degeneracy for very 
young and very old galaxies, allowing to statistically discrim¬ 
inate between the SMC and the Calzetti curve. 

According to our statistical test, the percentage of galaxies 
best fit with the SMC curve is < 10% of the MIPS sample 
(for any SFH), while galaxies best fit with the Calzetti curve 
represent ^25%. The rest of the sample does not have relative 
probability passing our threshold, and it represents ^ 65% of 


4.3. General properties of the sample 

In this section, we compare the properties of the MIPS sam¬ 
ple to those of the entire sample and with the properties of 
z ^ 2 star-forming galaxies taken from the literature in order 
to highlight the relevance of our method and the representa¬ 
tiveness of our sample. 

The stellar masses and SFRS derived assuming a constant 
SFH are consistent with the relation between M^ and SFR 


derived in iDaddi et al.l (l2QQ7h a ndiReddv et al.l (l2012bl) . with 
a rms scatter of 0.37 dex as in IReddv et al.l (12012bh (Figure 
El). The rms scatter for the rising and the declining SFHs 
are 0.56 and 0.75 dex, respectively. Using a Spearman cor¬ 
relation test, we find that SFR and M^ are weakly correlated 
(the strongest correlation is found for the declining SFH, with 
p = 0.22 and a = 2.56). We perform linear fits to our sample 
and subsamples (total MIPS sample and the MIPS-detected 
sample). While the stellar masses and SFRs for the MIPS 
detected sample do not depend strongly on the assumed star 
formation history, for galaxies without MIPS data the effect of 
star-formation history is more important. While ages derived 
with the rising star formation histories are similar to those de¬ 
rived assuming a constant SFH (A(logAge) = 0.0, a = 0.2), 
dust extinctions are higher relative to the ones derived with 
the constant SFH (A(logAv) = 0.2, a = 0.2). Accordingly 
the SFRs are also higher (A(logSFR) = 0.3, a = 0.3, see also 
ISchaerer & Pell6l[200l: Ide Barros et al.ll2014l) . For objects in 
the MIPS sample, the differences in SED parameters obtained 
with different SFHs are reduced relative to those derived when 
no IR constraints are available. 

While deriving the relation between SFR and M^ at z ^ 2 
is not the main focus of this paper, our method used to de¬ 
rive physical properties of galaxies leads to results consistent 
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0.5 


- 0.5 


- 16.5 - 16.0 - 15.5 

Log (ffiL/erg cm'^) [Observed] 



Figure?. Comparison between logfEL(obs) and log(fF jrohsi/fFJi S EPi) fo r 
the entire sample. Emission line fluxes are taken from lErb et al.l (I2006bn . 
IKulas et al.l and IKulas et all GoH), hereafter E06, K12, and K13 re¬ 

spectively. The blue line shows the median log(fEL(obs)/fEL(SED)) value and 
the grey area the la dispersion. 

Errors are determined by combining measurement error of the observed 
fluxes with the errors derived from parameter degeneracy (black cross in the 
right top corner). 

with what we expect for typical z^2 star-forming galaxies in 
terms of their SFRs and stellar masses, if we assume a con¬ 
stant SFH. 

To assess if our sample is biased in terms of dust proper- 
ties, we conipare ob served UV p slopes from the sample from 
iReddy et'd] (1201 5l) with the slopes of our sample in Figure [6l 
We derive the slopes using the filters probing the rest-frame 
range A = 1500-2500A (^e.g.. ICastellano et ^l2012h . Both 
our entire and MIPS s amples exhibit distributions similar to 
the iReddy et^ (1201 5h sample, with a small galaxy fraction 
probing redder (3 slopes (/3 > -0.5). 

Thus, our sample is likely representative of > IO^Mq z ^ 2 
galaxies in general. 

4.4. Reproducing emission line fluxes 

We now compare the emission line fiuxes predicted 
by our SED fitting code (see Section [3]) with the ob¬ 
served emission line fiuxes. We mainly study the re¬ 
sults for the Ha line because we have Ha measure¬ 
ments for 137 galaxies allowing for significant statis¬ 
tics. Furthermore, the [N ii] A6583 to 11(3 flux ratio 
depends strongly on metallicity with A6583/-^H/3 = 

0.175 for Z = 0.004 and ^[N ii] a 6583 /^h^ “ 0.404 for 
Z = 0.02, while the [O ill] A5007 to H/3 flux ratio is 
less metallicity dependent (^[Q iii] a5007/^h/ 3 = 4.752 and 
FjQ n il A 5 nn 7 /^H /3 = 4.081 for Z = 0.004 and Z = 0.02, respec¬ 
tively, (AnderT&TrhiewrAiyenslebe^llOOll), but more sen¬ 
sitive to the H II physical conditions li ke ionization parame¬ 
ter or leakage of ionizing photons (^e.g.. lKewlev et al.ll2013al: 
iNakaiima & Quchill2014)) . 

To quantify this comparison, we computed Afpp = 
log(fEL(obs)/fEL(SED)), where fEL(obs) is the observed emission 
line fiux and fEL(SED) is the SED predicted emission line fiux. 
We remind the reader that the SED fitting assumes that the 
color excesses of the nebular regions and the stellar contin¬ 
uum are equal, i.e., E(B - V)steiiar = E(B - V)nebuiar (derived as¬ 
suming the Calzetti extinction curve). Under this assumption, 
the Ha extinction is Ah^ = 3.33 x E(B - V)steiiar- The compar¬ 
ison between our predicted and observed fiuxes for the entire 


Figure 8. A log Iel vs. E(B-V)steiiar for the entire sample. The dashed line 
and the dotted line show the expected relation from equation [5] and equation 
[6| respectively. 



E(B-V)stellar 

Figure 9. Same as Fig.[8]for the MIPS detected sample. We show 2 different 
fits to the Ha data: in blue, we show a linear fit to the data, with the slope 
and intercept as free parameters. Additionally, in red, we show the linear fit 
where we fix the intercept to be zero. 


sample is shown in Figure [71 We assume a constant SFH be¬ 
cause this SFH provides the best match to the z ^ 2 galaxy 
physical properties (Sect. lOi . On average, our method used 
to model nebular emission from the photometry provides Ha 
fiuxes consistent with observations: the median shows a slight 
underprediction of the fiux by ^ 3% with a = 0.31 dex. 

We first study how Alogfpp correlates with the stellar color 
excess (Figured]). We compare our results with the relation- 
ship between nebuly and color excess for local galaxie^ from 
ICalzetti et al.l (120001) : 

E(B - V)stellar = (0.44 ± 0.03) X E(B - V)nebular (4) 

The extinction at a wavelength A is related to the color ex¬ 
cess E(B-V) and to reddening curve k(A) by Ax = k(X) x 
E(B-V). At the Ha wavelength (6562.SA), we have 
^Calzetti (Ha) = 3.33 and /:cardeiii(Ha) = 2.52, which leads re- 

^ Equation |4| is often misinterpreted as both color excesses being derived 
with the same attenuation curve, but the intent was that the nebular regions 
are described by a line-of-sight (e.g.. Milky Way) extinction curve (e.g., 
ISteidel et al.l[20ll : IReddv etaLll^^ . 
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spectively in terms of Ha extinction (assuming Equation [4]): 

AHa,Cardelli = 5.12 X E(B - (5) 

AHa,Calzetti = 7.55 X E(B ~ (6) 



• (2) We fit the MIPS detected sample assuming Z? = 0 in 
equation [T0| (red). 



Log (M*/Mo) 


Figure 11. A log Iel vs. M* for the entire sample. Typical error bar is shown 
in the top right corner. 


Figure 10. A log Iel vs. SFR for the MIPS detected sample. Linear fit to 
the data is shown with the blue line. 

Our comparison between predicted fluxes (based on 
SFRsed derived assuming energy balance) and observed 
fluxes should allow us to discriminate between equations [5] 
and[6l If the relative color excess of the UV continuum and 
the nebular emission lines at z ^ 2 is similar to the one derived 
in local starburst galaxies, we expect to overpredict emis¬ 
sion line fiuxes because we assume Ay, nebular = Ay, stellar- Fuf- 
thermore, if the nebular color excess is related to UV color 
excess as in local starburst galaxies (i.e., E(B - V)nebuiar = 
constant x E(B-V)steiiar)’ we should find an anticorrelation 
between A log fsL and E(B - V)steiiar- If we assume, as we do 
in our SED fitting procedure, that the same attenuation curve 
applies to both nebular and stellar emission, we get: 



fEL(obs) 

fintrinsic X 10-0-4x^A,obs, 

(7) 


fEL(SED) 

fintrinsic X 10^'4 x^a<sed, 

log 

f fEL(obs) \ 

V fEL(SED) J 

= -0.4 X (AA(obs) - Aa(SED)) 

(8) 

AlogfEL = 

= -0.4 X k(X) X (E(B - VU - E(B - 

(9) 


A Spearman correlation test performed on data from Figure [8] 
leads to p = -0.31 with standard deviation from null hypoth¬ 
esis cr = 3.63 (Table [B, which indicates that A log fpL and 
E(B-V)steiiar are possibly anticorrelated. We also compare 
the results expected from equation [5] and equation 0 with our 
data. While our method is able to reproduce observed emis¬ 
sion lines within a factor < 2, uncertainties remain too large 
to derive a possible relation between E(B-V)steiiar and A log 
fEL- To decrease the uncertainties, we focus on the MIPS de¬ 
tected sample (Fig. Hi. We compare the results expected from 
equation [5] and equation^with linear fits to our data: 

A log fEL = X E(B - V)stenar + ^ (10) 

and we perform two linear fits: 

• (1) We fit the MIPS detected sample (blue in Figure [8]) • 


Assuming b = 0 means that we consider that our SED fitting 
code is able to reproduce the emission line fiuxes for dust free 
galaxies. The intercept of fit (1) means that we underestimate 
observed emission line fiuxes for dust free galaxies by a factor 
- 2 . When assuming b = 0,E(B - V)stenar and E(B - V)nebuiar 
are similar (assuming a Calzetti extinction curve for both 
emission). This resu lt is consistent with the conclusion from 
IShivaei et al.l (l2Q15h . However, fit (1) shows that we are not 
predicting exactly emission fiuxes for galaxies with no extinc¬ 
tion. If we assume that there is a systematic offset between 
predicted and observed fiuxes, not dependent on E(B-V)steiiar, 
we can compare the slope of our relation with expectation 
from equation [5] and [6l Then, our slope is consistent with 
equation [21 

4.5. Relation between nebular color excess and other 
physical properties 

In this section we determine wether A log Iel possibly cor¬ 
relates with other physical properties such as SFR, M^, and 
sSFR. We summarize the correlation test and their associated 
significances for the different samples in Tabled] 

We find an anticorrelation between log SFR and A log fEL 
with a level of confidence > 6cr (p = -0.53, a = 6.67). This 
relation could explain the discrepancies among different pre¬ 
vious studies. We derive the relation between A log fEL and 
log SFR using the MIPS detected sample because uncertain¬ 
ties are too large when IR information is missing. Assuming 
a constant SFH, the relation from FigurefTOl implies: 

E(B - V)stenar “ F(B - V)nebular = 
(-0.50±0.10) X logSFR + (0.91 ±0.18) (11) 

This anticorrelation between A log fEL and log SFR implies 
that for z ^ 2 star-forming galaxies with log SFR < 
we have E(B - V)stebar ^ F(B - V)nebuiar- While uncertainty is 
large, this result is consistent with the SFR value derived in 
iReddv et^ (l2015h . We discuss this point further in Section 
(531 

We find no correlation betwen M^ and A log fEL IFig.fTTI). 
with p = 0.16 and a = 1.86. On the contrary, we find an an¬ 
ticorrelation between sSFR and A log fEL with p = -0.52 and 
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(7 = 6.48 (Fig.IEli. Again, the uncertainties are too large to 
derive a relation between A log fsL and log sSFR for the en¬ 
tire sample, and we use the MIPS detected sample to derive a 
relation. From Figure [121 we get: 


E(B - V)steiiar “ ^(B - V)nebular “ 

(-0.31 ±0.07) X log(sSFR)-r (0.16 ±0.03) (12) 



Figure 12. A log Iel vs. sSFR for the MIPS detected sample. Linear fit to 
the data is shown with the blue line. 


We conclude that the analysis of our sample sup¬ 
ports marginally the relation derived between E(B - V)steiiar 
and E(B-V) Eiiior in local star-forming galaxies (Eq. |5l 
ICalzetti et al.ll200Ql) . While we do not find a correlation be¬ 
tween the stellar mass and E(B-V)steiiar“E(B-V)nebuiar’ we 
find anticorrel ations between SFR /sSFR and E(B - V)stebar “ 
E(B-VU„i,, (lReddvetal.1120151) . 

5. DISCUSSION 

5.1. Modeling of nebular emission 

One result of our study is that our SED fitting code is able to 
reproduce consistently the broad-band photometry and emis¬ 
sion line fluxes, within their respective uncertainties, if we as¬ 
sume energy balance. This result is not a trivial finding as our 
modeling of the nebular emission relies on several assump¬ 
tions regarding the Lyman continuum escape fraction, ISM 
physical conditions, empirical ratios between emission lines, 
and the IME. 

However, we can not rule out that a different set of phys¬ 
ical properties could produce similar emission line fluxes. 
Eor example, the Ha flux woul d decrease by a factor 1.8 
if we consider a Chabrier IME (IChabried l2003h instead of 
a Salpeter IME. Eurthermore, the absolute fiux of hydro¬ 
gen emission lines are directly proportional to the number 
of ionizing photons per sec ond, and so to the L yman con¬ 
tinuum escape fraction te.g.. lKrueger et al.lll995l) . Neverthe¬ 
less, we show in Eigure[T3]the comparison between predicted 
and observed [O ill] A5007 and [N ii] A6583 fluxes. We 
have shown in the previous sections that the effect of dif¬ 
ference between nebular and stellar color excesses is neg¬ 
ligible in the range of SER (and sSER) probed with our 
sample. Therefore, A log fpL for these two lines is mostly 
sensitive to the metallicity for [N ii] A6583 (from Zq to 
O.2Z0 the [N II] A6583/Ha ratio decreases by a factor 2.34, 


lAnders & Eritze-v. Alvenslebenll2003l) . and to the ISM phys - 
ical conditions for [O III] A5QQ7fe.g.. iKewlev et al.ll20T^ . 
Our method is able to reproduce these lines with a similar ac¬ 
curacy as it reproduces Ha, with a median A log fsL^ 0 and 
a = 0.34 dex. 



Figure 13. Same as Figure [2]for the [N ii] A6583 and [O iii] A5007 mea¬ 
surements. The blue line shows the median A log Iel value and the grey area 
the la dispersion. 



O Ha E06 
• HaK13 

7 8 9 

Log (r^/yr) 


Figure 14. A log Iee vs. timescale rj for the MIPS detected sample assum¬ 
ing a declining star-formation history. We show the result only for the Ha 
data for more clarity. 


We show that the method used in this work to model 
nebular emission (ISchaerer & de Bari^l2009L 120101) is able 
to reproduce Ha, [N II] A6583, and [O III] A5007. This 
is an important result because of the implications for 
higher redshift studies (z > 3), where emission lines can 
not be measured with current instruments while observed 
SEDs exhibit evidences of impact of nebular emission (e.g. . 
Stark et aP 120131: Ide Barros et al .11201 dUNavyeri et al.l 12014 

Marmol-Queralto et~^l2015l: ISmit et al.ll2015h. 

5.2. Star-formation history 

Several studies have used the comparison between direct 
star-formation tracers and SER(SED) to put constraints on 
star-formation timescales. Their main conclusion is that de¬ 
clining SEHs result in SERs(SED) that are inconsistent with 
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Table 1 

Spearman correlation test between A log fsL and physical parameters 


Physical parameters MIPS detected (A/^ = 14) MIPS sample (A/^ = 34) Entire sample (A/^ = 137) 

p a p a 


E(B V) stellar 

-0.46 

1.61 

-0.35 

1.96 

-0.31 

3.63 

SFR 

-0.62 

2.31 

-0.78 

5.68 

-0.53 

6.67 

M* 

0.66 

2.58 

0.13 

0.73 

0.16 

1.86 

sSFR 

-0.76 

3.21 

-0.57 

3.49 

-0.52 

6.48 


^ Spearman rank correlation coefficient. 

^ Standard deviations from null hypothesis. 


those inferred from multi-wavelength tracers of star forma¬ 
tion (e.g., the sum of the IR and UV), or that it is necessary 
to assume a lower limit on the decay timescale (r^ > 300 Myr; 
IWuvts et alJlMItiReddv etani2012bl:fPrice et al.ll20Tl) . To 

gain insight on the star-formation timescales for both rising 
and declining SFHs, we rely on the MIPS detected sample 
since the SFR uncertainties are much lower when L(IR) is 
used to perform SED fitting assuming energy balance (the un¬ 
certainties derived from SED fitting are 0.1 dex versus 0.35 
dex without assuming energy balance). Therefore, we use 
both IR a nd emission lines to place constraints on the SFH 
timescale (ISchaerer et al.ll201^ . 



Figure 15. A log fsL vs. the ratio between age and timescale r/r^ for the 
MIPS detected sample assuming a declining star-formation history and and a 
relation between E(B=V)nebuiar and E(B=V)steiiar described by Equation [TT] 
The dash-dotted line shows r/r^ = 1. We show the result only for the Ha data 
for more clarity 


Since we assume energy balance, the IR luminosity pre¬ 
dicted from SED fitting are consistent with L(IR) 24 ^m by con¬ 
struction. The difference between predicted and observed 
emission lines fiuxes assuming a declining SFH do not dif¬ 
fer significantly from other SFHs, except for galaxies with 
blue stellar color excess. In Figure [141 we show the re¬ 
lation between A log fpL and the timescale rj and in Fig¬ 
ure \T5\ the relation between A log fpL and the ratio between 
age and timescale f/r^. Overall, the SED fitting reproduces 
the observed Ha fiuxes with the same accuracy as for the 
entire sample (i.e., the fiuxes are reproduced within a factor 
2). There are two galaxies for which the fiuxes are signif¬ 
icantly underpredicted (A log fEL>0.5) and they are fit with 
short timescales (< lOOMyr). In Figure [151 we see that the 
age to timescale ratio can be as high as ~ 5 with the pre¬ 
dicted emission line fiuxes still consistent with the observed 
values. The two galaxies for which we significantly underpre¬ 


dict the emission line fiuxes have the lowest IR luminosities 
of the sample (< IO^^Lq). If we add the Ha fiux measure¬ 
ment constraints to the SED fit for these two galaxies (ex¬ 
cluding solutions providing Ha fiuxes inconsistent with the 
observed ones within la), timescales are > 500Myr and 
the age to timescale ratios are f /r^ < 1. Our conclusions re- 
garding declining SFHs are consiste nt with the literature (e.g., 
IWuvts et n]l201 ll: iPrice et aO 120141) but for a fraction of our 
MIPS detected sample (~ 1/3), we find that galaxy SEDs are 
fit with short timescales (30Myr < rj < IGyr) and relatively 
large tjrd ratios (0 < < 5), and have predicted physical 

quantities (L(IR) and Ha fiux) consistent with observations. 

We overpredict the emission line fiuxes for objects fit with 
a short exponential timescale (r^ < lO^yr) for the rising SFH 
(data not shown), while we underpredict fiuxes for fit with a 
large timescale (> lO^yr). 



Log (SFRh./M^ yri) 


Fi gure 16. SERnr. y s. S FRuy assuming E(B-y),,^ii, = E(B-V)nphniar for 
the lErb et'dl J2006ar) and IForster Schreiber et al.U200^ samples. The hori¬ 
zontal red line sh ows log SFR =1.82 (Se ction l4.5l EQuationfTll. the blue line 
lo^ SFR - 2 25 jy)shikawa et al.ll201Cir) . and the green line log SFR =1.56 
iReddv et al .11 20 151 converting the value from a Chabrier to a Salpeter IMF). 
The dash-dotted line is the expected log SFR(Ha)-log SFR(UV) relation as¬ 
suming that E(B - V)nebuiar “E(B - V)stenar related to log SFR as described 
in Equation[TT] The dashed line is the one to one relation. 


To perform SED fitting and get results consistent with ob¬ 
served emission line fiuxes, assuming energy balance (and so 
having IR information), we show that the timescale used with 
a declining SFH is rj > 30Myr and that the age to timescale 
ratio is < 5. For a rising SFH, the timescale is in the range 
300Myr < < IGyr. 

5.3. Comparison with other studies 

We find an anticorrelation between SFR and A log fpL 
which we interpret as a correlation between SFR and 
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E(B - V)steiiar “E(B - V)nebuiar (Equation[TT]): the nebular color 
excess becomes larger than the s tellar color excess with in¬ 
creasing SER. As pointed out by iReddv et^ (I201QI) . previ¬ 
ous studies included galaxies with different ranges in SER 
and M^. Therefore we check now if our relation found be¬ 
tween SER and A log fsL could explain the discrepant results 
on this matter. 

We compare the lErb et al.l (l2006ah sample, for which 
the nebular and stellar color excesses are similar, and the 
lEorster Schreiber et all (l2009l) sample, for which the nebu¬ 
lar color excess is higher than the stellar color excess (Equa¬ 
tion H]). It is worth noting that a fraction of our sample is 
drawn from the lErb et"^ (l2006ah sample. 

We show in Eigure [16] the comparison between SERuv 
and SERhq. for these two studies assuming E(B - = 

E(B-V)nebuiar’ ^ Calzetti attenuation curve, and after con¬ 
verting the SERs from a Chabrier IME to a Salpeter IME. 
We expect to see an increasing discrepancy between SERhck 
and SERuv above the value derived from Equation [TT] with 
lo g SER = 1.82. We also show similar SER thresh olds found 
in lYoshikawa et all (1201 Ql) and (iReddy et al.ll2015h . 

Because we do not have Baimer decrement measure¬ 
ments in order to correct Ha, we consider the dust-corrected 
SERuv as the one with which we compare SERua- We 
can predict the difference between SER(Ha) derived assum¬ 
ing E(B-V),teiiar = E(B-V)nebuiar and SER(Ha) derived as¬ 
suming the relation between SER and the difference in the 
color excesses as in Equation [TT] (dash-dotted line in Eig¬ 
ure [Ell. At log SER= 1.82, this difference is zero and 
reach 0.56 dex at SER(UV)con-ected = lO^M 0 yr“^ This ex¬ 
plains why we do not see any significant discrepancy between 
SER(U V) and SER(Ha) i n Figur e [T^ below the SER limit de¬ 
rived in lYoshikawa et al.l (l 2010 l) . Assuming E(B - V)steiur = 
E(B - V)nebuiar l^ads to underpredict SFRu n. in comparison 
with S ERuv above the threshold found in lYoshikawa et al.l 
(I 2 OIOI) . Therefore we conclude that the discrepant result 
about the difference between nebula r and stellar color ex¬ 
cesses is likely a selection effect, the iForster Schreiber et al.l 
(I 2 OO 9 I) sample probing higher SERs. The apparent discrep¬ 
ancy b etween Equation [U] and the iForster Schreiber et al.l 
(I 2 OO 9 I) sample can be explained b y the large dispers ion in 
color excess for individual galaxies (iReddv et al.ll2015h . 

6. CONCLUSIONS 

We use a sample of 149 star-forming galaxies spectroscop¬ 
ically confirmed at z ^ 2 with optical emission line measure¬ 
ments to investigate the relative dust attenuation of the UV 
continua and the nebular emission in these galaxies. For a 
subsample of 41 galaxies, MIPS 24/im data are available and 
we use them to derive the infrared luminosity of these galax¬ 
ies. We compare all the available observed/derived proper¬ 
ties of our sample (emission line fiux, IR luminosity) with 
predictions from stellar population synthesis models, taking 
into account nebular emission, and assuming energy balance 
between UV and IR emission. We explore a large param¬ 
eter space, exploring different star-formation histories (con¬ 
stant, declining, rising) and the impact of different attenuation 
curves (Calzetti and SMC curves) on our results. 

We compare the properties (stellar mas ses, SERs) of our 
sample with proper ties found in literature (iDaddi et al.ll2007l: 
iReddy et^ l2012ah to ensure the representativeness of our 
sample. It appears that we are not able to reproduce the ob¬ 
served UV p slope for the ULIRGs in our sample when we 
assume energy balance and we exclude these 7 galaxies from 


our analysis. Our results confirm that the Calzetti curve is 
not appropriate for ULIRGs (Appendix and an appropri¬ 
ate attenuation curv e should be define t o derive properly their 
physical properties (IReddv et al.ll201Ql:[Casev et al.ll2014l) . 

Our conclusions are as follow: 

• Assuming energy balance, our SED fitting code is 
able to reproduce the observed emission line fluxes for 
galaxies in our sample (within a factor < 2). 

• Young galaxies (< 100 Myr) are statistically better 
fit with an SMC curve rather than a Calzetti curve 
(IReddv et 

• We find a correlation between SER and 
E(B-V).tpiiar - E(B-V)_u„u^, and also between 
sSFR and E(B-V),teiiar " E(B-V),ebuiar. The color 
excess difference increases with increasing SER and 
sSFR. The correlation between the amount of extra 
attenuation toward nebular emission and the SER could 
explain discrepancies among previous studies about 
the n ebular attenuation relative to stellar attenuation 
fe.g.. lErb et al.l l2006al: IForster Schreiber etd] 120091: 
IReddv et al.ll2'(m 

• We use SED, IR and emission line fiux measurements 
to constrain the star formation history. Overall, a large 
range of SFHs can reproduce the observables. We find 
that assuming a declining SFH, the timescale is rj > 
30Myr and the age to timescale ratio is t/rd <5. For a 
rising SFH, the timescale is 300Myr <Tr< IGyr. 

We look forward to comparing our re sults with cur¬ 
rent ongoing survey (IReddv et al.l l2015h . This will 
also allow to determine much more precisely the ratio 
E(B - V)stenar/E(B - V)nebuiar’ determine if there is an evo¬ 
lution of this ratio from low to high redshift. However, our 
results suggest that this kind of observations must aim a large 
range of SERs and stellar masses to avoid any bias. It would 
also be interesting to measure B aimer decrement of ULIRGs 
to derive the best method to infer their physical properties 
such as stellar mass, SER, and dust attenuation. 
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for suggestions that clarify the text and the analysis. We ac¬ 
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the data presented here. We thank Alice Shapley, Valentino 
Gonzalez, Pascal Oesch, Bahram Mobasher, Daniel Schaerer, 
Panos Sklias, Alberto Dominguez, and Matthew Hayes for 
useful discussions about this work. We also thank Roser 
Pello, Daniel Schaerer, and Panos Sklias for their support with 
HyperZ. 

A. ULIRG PHYSICAL PROPERTIES 

As shown in Figure [2] for the ULIRGs, the SERs derived 
from SED fitting are inconsistent with those derived from UV 
and IR luminosity. The mismatch between SFRuv+ir and 
SFRsed arises from the inability of our SED fitting code to 
reproduce both the derived IR luminosity and the observed 
UV (3 slope. UV (3 slope can be used to derive the dust at¬ 
tenuation assuming an attenua tion curve for typic al (i.e., L*) 
galaxies at high-redshift (^e.g.. IReddy et al.ll2QlQl) . We show 
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Figure 17. Observed SEDs (black thick points) of the seven ULIRGs in our sample with the best-fits for three models: constant SFH, age > 50Myr, and 
Z = [0.004,0.02] (blue); constant SFH, age > 50Myr, Z = [0.004,0.02], and assuming energy balance between UV and IR (red); and declining SFH, age free, 
Z = [0.0004,0.004,0.02], and assuming energy balance (black). Crosses show the synthesised flux in the Alters. For each model and galaxy, we show the the 
UV /3 slope associated with the best-fit solution, and the predicted infrared luminosity when relevant. 


in Figure [T7] the SED fits of the seven ULIRGs in our sam¬ 
ple for three different models: two assuming a constant SFH, 
with and without assuming energy balance, and one assuming 
a declining SFH with relaxed assumptions on age and metal- 
licity, and assuming energy balance. This latter model is used 
because it is the model with the maximum number of degrees 
of freedom and should be more able to reproduce the SED. 

The observed P slopes are bluer that what we would ex- 
pect given the high IR luminosities, assu ming a Calzetti curve 
(iMeurer et al .111 9991 : ICalzetti et al.ll200Ql) . This is also consis¬ 
tent with the better agreement between SFRsed and SFRuv+ir 
for the ULIRGs assuming a declining SFH and a large range 
of ages and metallicities: the p slope is more sensitive to 
change in age an d metallicity for a declining SFH than for 
a constant SFH (iLeitherer & HeckmanI Il995() . Indeed, re¬ 
laxing assumptions on age and metallicity for the declining 
SFH leads to slightly younger ages (age > 30Myr) and SED 
fits with the minimum allowed metallicity (Z=0.0004), which 
lead to better reproduce the observed P slope, but with metal¬ 
licity values unrealistically low (e.g., ?). However, even this 
last set of assumptions does not yield satisfactory fits of the 
UV p slope for the most IR luminous galaxies: the P slopes 
predicted by the SED fits are redder than the observed ones. 

Our study confirms previous results regarding ULIRGs 


(e.g., iGoldader et all 1200^ iReddv et al.l l2^ : ULIRGs 
generally have signficantly redder IR/UV luminosity ratios 
that would be infer r ed from their UV s lopes assuming the 
iMeurer ef^] (1 19991) : ICalzetti et an(l200Ql) attenuation curves. 

Furthermore, depending on whether we assume energy bal¬ 
ance, we have AlogfEL = log(fEL(obs)/fEL(SED)) as high as 
^1.0 dex. This is due to the large difference in dust attenua¬ 
tion derivation found for the ULIRGs between the two method 
used here (with or without energy balance): the mean differ¬ 
ence is A Ay = Ay, energy balance ~ Ay^ no energy balance “ L6. 
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Figure 17. Continued. 
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